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An extensive study of the effects of correlations on both charged and neu- 
tral current weak interaction rates in dense matter is performed. Both strong 
and electromagnetic correlations are considered. The propagation of particle- 
hole interactions in the medium plays an important role in determining the 
neutrino mean free paths. The effects due to Pauli-Blocking and density, spin, 
and isospin correlations in the medium significantly reduce the neutrino cross 
sections. Due to the lack of experimental information at high density, these 
correlations are necessarily model dependent. For example, spin correlations 
in nonrelativistic models are found to lead to larger suppressions of neutrino 
cross sections compared to those of relativistic models. This is due to the 
tendency of the nonrelativistic models to develop spin instabilities. Notwith- 
standing the above caveats, and the differences between nonrelativistic and 
relativistic approaches such as the spin- and isospin-dependent interactions 
and the nucleon effective masses, suppressions of order 2-3, relative to the 
case in which correlations are ignored, are obtained. Neutrino interactions in 
dense matter are especially important for supernova and early neutron star 
evolution calculations. The effects of correlations for protoneutron star evolu- 
tion are calculated. Large effects on the internal thermodynamic properties of 
protoneutron stars, such as the temperature, are found. These translate into 
significant early enhancements in the emitted neutrino energies and fluxes, 
especially after a few seconds. At late times, beyond about 10 seconds, the 
emitted neutrino fluxes decrease more rapidly compared to simulations with- 
out the effects of correlations, due to the more rapid onset of neutrino trans- 
parency in the protoneutron star. 

PACS numbers(s): 13.15.+g, 26.60.+C, 97.60.Jd 
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I. INTRODUCTION 



The calculation of neutrino interactions in hot, dense matter is highly relevant to the 
study of supernovae and protoneutron stars. Neutrinos are thought to be intimately involved 
in the supernova explosion mechanism. In addition, the timescales over which protoneutron 
stars deleptonize and cool, which are crucial in predicting the neutrino light curve from a 
supernova, are determined by neutrino opacities. In an earlier paper [ij, we developed a 
formalism to determine both neutral and charged current opacities for interacting matter, 
accounting for in-medium mass and energy shifts given by its underlying equation of state 
(EOS). The effects of degeneracy and relativity were incorporated, and nuclear interactions 
modelled by both nonrelativistic potential and relativistic field-theoretical approaches were 
studied. In this paper, we extend this treatment to include important sources of suppression 
and enhancement due to in-medium correlations. In particular, we include effects due to 
both strong and electromagnetic correlations in the system. While neutrinos can couple 
to many-particle states and induce multi-pair excitations, both of which could modiy the 
response function developed in Ref. |L], single-pair excitations dominate over multi-pair 
excitations for the kinematics of interest in neutrino scattering and absorption. Thus, we 
focus on the effects of correlations on the single pair excitation spectrum. 

Pioneering work in calculating neutrino mean free path in uniform nuclear matter was 
performed by Sawyer who showed that the effects due to strong interactions are important 
and that the relation between the EOS and long wavelength excitations of the system may be 
exploited to compute the weak interaction rates consistent with the underlying dense matter 
model. Iwamoto and Pethick || computed the neutrino mean free paths in pure neutron 
matter for the densities of relevance for neutron stars by calculating the dynamical response 
within the framework of Fermi-liquid theory • Horowitz and Wehrberger || studied the 
influence of correlations on the neutrino scattering in dense matter in a relativistic field- 
theoretical model. More recently, Fabbri and Matera || have calculated neutrino scattering 
in asymmetric nuclear matter including the effects of the exchange interaction. 

These earlier works on neutrino interactions in dense matter have shaped our approach. 
First accounts of our results may be found in Refs. []ID|,|rT| . We urge the reader to consult 
the simultaneous and independent work of Burrows and Sawyer in Ref. |12[ , with which our 
work has some overlap. However, important differences, both in the microphysical inputs to 
the opacity calculations and in the underlying EOS models, distinguish our approach and 
results from theirs. We will present a unified approach to calculate both the neutral current 
scattering and charged current absorption reactions for the composition and temperatures of 
relevance to the PNS evolution. A systematic analysis of the various nuclear physics inputs 
such as the particle-hole (abbreviated p-h) interactions in density, spin and isospin channels 
is undertaken keeping in view the known nuclear ground state and excited state properties. 
We will isolate the important sources of enhancement and/or suppression due to in- medium 
correlations and highlight the important role of spin and isospin dependent p-h interactions. 
The need for further efforts to pin down the poorly known spin-dependent interactions is 
clearly brought out by our investigations. 

The neutrino cross sections and the EOS are intimately related. This relationship is 
most transparent in the long- wavelength or static limit, in which the response of a system 
to a weak external probe is completely determined by the ground state thermodynamics 
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(EOS). Thus, in this limit, neutrino opacities which are consistent with a given EOS || 
can be calculated. However, when the energy and momentum transferred by the neutrinos 
are large, full consistency is more difficult to achieve. Particle-hole and particle-particle 
interactions [pT0|— [121 , multi-pair excitations [13], and the renormalization of the axial charge 



I1] , |T5|| are examples of effects that may influence both the EOS and the opacity. In this 



work, we investigate the influence of particle-hole excitations utilizing the Random Phase 
Approximation (RPA), in which ring diagrams are summed to all orders. This may be viewed 
as the first and minimal step towards achieving self-consistency between the response of 
matter to neutrino-induced perturbations and the underlying EOS, similar to what has long 
been established in the study of collective excitations in nuclei (multipole giant resonances) 
through electromagnetic probes Jl6| . 

Throughout this paper, we employ both nonrelativistic potential and relativistic field- 
theoretical models for the EOS and the evaluation of the neutrino-matter cross sections. 
This is done in order to evaluate effects that are generic as well as effects that are due to the 
precise description of the nucleon-nucleon interactions. In Sec. II, we study the response 
of pure neutron matter and calculate the neutrino cross sections to gauge the effects due 
to RPA correlations. In Sec. Ill, we calculate the influence of many-body correlations on 
the scattering and absorption reactions in both nuclear matter and mult i- component stellar 
matter (charge neutral, /5-equilibrated matter). In Sec. IV, we compute protoneutron star 
evolutions in order to evaluate the effects of correlations. Our conclusions and outlook are 
contained in Sec. V. Explicit formulae for the various p-h interactions and polarization 
functions needed to calculate the neutrino cross sections are collected in Appendices A and 
B, respectively. 



II. NEUTRAL CURRENT NEUTRINO CROSS SECTIONS IN NEUTRON 

MATTER 

We begin with pure neutron matter to illustrate the various strong interaction correlation 
effects that influence the neutrino mean free paths in dense matter. 

A. Nonrelativistic Models 

In the nonrelativistic limit for the neutrons, the cross section per unit volume or the 
inverse mean free path for the neutral current reaction u + n — > v + n is given by || 

El [1 -/,(*)] 



V d£l 2 dq 8tt 3 



x 



-v 



[l + cos6)Soo(q ,q) + c 2 A (3-cos9)S w (q ,q) , (1) 



where cy = 0.5 and ca = 0.615 are the neutral current vector and axial vector couplings of 
the neutron and G 2 F = 6.89 x 10 -10 MeV _5 m _1 , Ei and £"3 are the incoming and outgoing 
neutrino energies, qo = Ei — E 3 and q = ki — k 3 are the energy and momentum transfers 
in the scattering reaction. The neutron matter response, which depends on the kinematical 
variables q and q = \q\, is characterized by the function S(q , q) and is called the dynamic 
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form factor. Since neutrinos couple to the vector and axial vector currents of the neutron, the 
response is given in terms of the density-density correlation function Soo(<70) q) an d the spin- 
density correlation function S w (qo, q). The subscripts appearing in the response function 
S ar represent the total spin and isospin transferred along the p-h channel and are labelled 
as follows: r = and a = for the neutral current density response, r = and a = 1 for 
the neutral current spin-density response, r = 1 and a = for the charged current density 
response, and r = 1 and o = 1 for the charged current spin-density response. In the absence 
of strong interaction correlations, S 00 = Si = S°(q ,q). The function S°(q ,q) is directly 
related to the lowest order polarization function IT (also known as the Lindhard function 
T7J ) through the principle of detailed balance, and accounts for correlations arising due to 



Pauli blocking llpf. Explicitly, 

Im n°(g , q) 



S°(qo,q) 



1 - exp (-q /T) 



f(E p )[l - f(E w+ql )) [1 - f(E p )]f(E w+ql ) 



q + E p - E\$+<ft + ie q + E p - E\p±$ - it 



(2) 



where f(E) = (1 + exp (E — yU n )/T) _1 is the Fermi-Dirac distribution function, /i n is the 
neutron chemical potential, and E p = p 2 /2M is the noninteracting nonrelativistic dispersion 
relation for the neutrons. The imaginary part of the polarization function may be explicitly 
evaluated when q < 2kp and qo < q, and is given by fllpOl 



Im U°{q ,q) 



M 2 T 



% ln ( 1 + exp[(e„ -/i 2 )/T] \ 
T n l v l + exp[(e_ + g -^4)/T] y ' 



2vrg 

P 2 - 1 (<7o - q 2 /2M) 



2M 4 q 2 /2M 



(3) 



The total cross section per unit volume or the inverse collision mean free path of neutrinos 
due to scattering is obtained by integrating over the qo — q space. 

The simplest modification to S(q , q) due to strong interactions arises due to the medium 
modification of the single particle dispersion relation In the mean field approximation 
(in which particles are assumed to move independently of each other in a common potential), 
the single particle spectrum may be cast into the form 

E(p) = p 2 /2M* + U, (4) 

in schematic models which are designed to reproduce the results of microscopic calculations. 
The single particle potential U and the (Landau) effective mass M* are in general density 
dependent. Because the functional dependence of the spectra on the momenta is similar 
to that of the noninteracting case, the only modification to the free gas results arises due 
to the nucleon effective mass. Thus, to obtain the mean field or Hartree response we need 
only replace the free nucleon mass M in Eq. (|3|) by the in-medium mass M*. The density 
dependence of M* is model dependent and in nonrelativistic approaches arises due to the 
presence of a momentum dependent interaction. Modifications due to M*{jib) are important 
and always act to decrease the scattering cross sections. This is chiefly due to the fact that 
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a lower M* decreases the density of states Nq = M*kp/Ti 2 near the Fermi surface and 
thereby the number of target neutrons available for scattering. Kinematically, a lower M* 
favors larger energy transfers by shifting the strength from low go to regions of high q Q . 
For thermal neutrinos, the dominant contribution to the total cross section comes from the 
region |go| < 7rT, since factors arising from the principle of detailed balance and final state 
blocking for the neutrinos exponentially suppress the phase space when |go| > tiT . Thus, 
for partially degenerate matter and for neutrino energies of order T, this shift in the strength 
leads to significant suppressions in the total cross sections. 

The results in Fig. [TJ show how M*(ub), which depends on the strong interaction model, 
modifies the neutrino mean free path as a function of density for different temperatures. For 
the illustrative results presented here, we have used M*/M = [H-a^ng/no)] -1 , characteristic 
of a large class of Skyrme models for tie < (3 — 4)no, where uq = 0.16 fm -3 is the nuclear 
saturation density. Denoting the value of the nucleon mass at n Q by Mq, the factor a = 
(M — Mq)/Mq. For the results shown in Fig. |T|, we have set Mq = 0.8. Our findings 
indicate that the neutrino mean free paths are roughly enhanced by a factor proportional to 
(M/M*) 2 . With increasing temperature, this enhancement decreases only moderately. We 
are thus led to the conclusion that for the densities and temperatures of interest, effects due 
to a dropping in-medium nucleon mass significantly increases the neutrino mean free path. 

We turn now to address the effects of strong interaction correlations. We first note that 
the nucleon-nucleon interaction for densities of relevance is not amenable to a perturba- 
tive treatment [pTp2f . Thus, we are forced to adopt approximate schemes to incorporate 



nonperturbative effects. The random phase approximation (RPA), also known as the ring 
approximation, is particularly suited to describe the response of matter at high densities. In 
RPA, ring diagrams are summed to all orders in order to incorporate p-h correlations. Sin- 
gle pair excitations created by the neutrinos are now allowed to propagate in the medium 
due to the presence of p-h interactions. The strength of the interaction in a specific p-h 
channel will therefore play a crucial role in determining the magnitude of the RPA correc- 
tions. Deferring to later sections the discussion of the interaction potentials, we set up here 
the basic formalism to calculate the cross sections for a given interaction Voo in the spin 
independent channel and Viq in the spin dependent channel. The RPA polarization, as a 
result of summing ring diagrams to all orders, is [18| 

n_° 
i 11,11° 



nf A = t^t™- (5) 



•i 3 

The dynamic form factor which incorporates these correlations is given by 



1 



Im n°(g , q) 



e i,j 



Sij(qo,q)= , , 

[1 - exp(g /T) 

Mffo, ?) = [!- VijRe U°(q , q)} 2 + [T^-Im n°(g , q)} 2 . (6) 

In addition to the imaginary part of the polarization, which is sufficient to evaluate the free 
gas response, the real part is also needed to evaluate the corrections due to strong interaction 
correlations. At finite temperatures, the real part is easily evaluated numerically by using 
the Kramers- Kronig relation 

Ren» ta , 9 ) = -ipr A, Imn ° ( ^> . (7) 

71 J-oo u) — q 
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The dielectric screening function e^j is the only modification that arises in RPA. If the 
interaction is repulsive €ij(q = 0,q) > 1, the medium response is suppressed, while an 
attractive interaction will result in an enhancement. If the repulsive interaction is strong, 
the RPA response predicts the existence of collective excitations such as zero-sound and spin 
zero-sound which could enhance the response when the energy transfer is finite and equal 
to the energy of the collective state. 



1. Particle-Hole Interactions 

The particle-hole (abbreviated as p-h) interaction, V%j, is in general a function of density, 
temperature, q , and q. For the temperatures of interest in the PNS evolution, q and q 
are typically smaller than a few tens of MeV's. Since the p-h interactions due to strong 
interactions are short ranged (~ 1/meson mass) compared with the typical wavelength of 
the excitations probed by the neutrinos, the momentum dependence of the interaction will 
not play a significant role. It is well known that the effective interaction of meson exchange 
models result in very large matrix elements due to the presence of a large short-ranged 



component and cannot be used in perturbative expansion [p2lp2|| . In nonrelativistic models, 
there are essentially two complementary methods that have been used to estimate the p-h 
interactions in nuclei and in bulk nuclear matter. We outline the basic features of these 
methods and employ the corresponding p-h interaction potentials to compute the neutrino 
scattering cross sections. 



Schematic Potential Models 



Skyrme-type models start from a zero-range, density, and momentum dependent po- 
tential from which an energy density functional is constructed. The force parameters are 
determined empirically by calculating the ground state in the Hartree-Fock approximation 
and by fitting the observed ground state properties of nuclei and nuclear matter. Using 
Landau Fermi-liquid theory, the effective p-h interaction is given by the second functional 
derivative of the total energy density with respect to the densities taken at the Hartree-Fock 
solution: V!y = (5 2 E(n)/SniSrij), and is usually expressed in terms of Fermi-liquid parame- 
ters. When the interaction is short ranged compared to the wavelength of the excitations of 
interest, justifying the assumption that only / = terms are retained, assumes a simple 
form. For example, in symmetric nuclear matter the p-h interaction is given by 

Vfa, k 2 ) = N Q l [F + G Q a x ■ a 2 + (Fq + G' a 1 - a 2 )n ■ r 2 ] , (8) 

where a and r are spin and isospin Pauli matrices, which act on the p-h states. To achieve full 
generality, Eq. (|8[) should be supplemented with contributions from tensor interactions. The 
Fermi-liquid parameters F ,Go,Fq, and G' are dimensionless numbers and N = 2M*/ci?/7r 2 
is the density of states at the Fermi surface. In neutron matter, since we have frozen the 
isospin degree of freedom, only the first two terms contribute (with Nq 1 = M*kp /it 2 ), and 
the effective interaction Voo = N^Fq and V w = N ~ 1 Go. The nucleon effective mass M*, 
which arises due to a momentum dependent interaction, is given by M* /M = (1 + iq/3). 
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Skyrme models have been successful in describing nuclei and their excited states ||16|| . In 
addition, various authors have explored its applicability to describe bulk matter at densities 
of relevance to neutron stars. The effective nucleon-nucleon interaction of the standard 
Skyrme model is given by the potential f2F 



V NN {r) = t (l + x P a )8{r) + \hn'{l + x 3 P a )5{r) 

o 

+ h 1 (l + x 1 P a )[k 2 6(r)+k^6(r)}+t 2 (l + x 2 P a )k^(r)k, (9) 

where, the force parameters t ,ti,t 2 ,t 3 , and the exchange force parameters x ,x 1 ,x 2 , and 
x 3 are empirically determined. The three-body interaction, which is written as a density 
dependent two-body interaction (through the parameters £3 and 7) and the momentum 
dependent interactions arising due to t\ and t 2 terms dominate the high density behavior of 
these models. 

The Fermi-liquid parameters for neutron matter, expressed explicitly in terms of the 
Skyrme parameters, are collected in Appendix A. Fig. ^| shows the results for four different 
sets of Skyrme parameters for Auq. A striking feature is the large density dependence of 
the Fermi-liquid parameters. While the qualitative behaviors of all four parameter sets are 



similar, significant quantitative differences exist. The models SkM* and SGII [^4j have been 
constrained by fitting the properties of systems with very small isospin asymmetries, while 
the models SLy4 and SLy5 were further constrained to reproduce the results of microscopic 



neutron matter calculations p5[ . All four models become unstable to spin oscillations and 
are driven towards a ferromagnetic ground state since Gq > —I for densities in the range 
(2 — 4)no- For > Go > — 1 > the RPA corrections enhance the spin response approximately 
by a factor (1 + G*o)~ 2 resulting in small neutrino mean free paths. For Go < — 1, the spin 
symmetric ground state is unstable. At the phase transition density, RPA response functions 
of the spin symmetric state diverge, and the formalism developed thus far is inapplicable. 
The response functions of the energetically favored, spin polarized ground state needs to be 
calculated. We do not attempt to do this here; we merely note that large enhancements in 
the RPA response functions are precursors to a phase transition. We present results only for 
riB < 2no, for which the ground state is spin symmetric. Spin instability is a common feature 



associated with a large class of Skyrme models [gfj j but is not realized in more microscopic 
calculations. It must be emphasized that the interaction in the spin dependent channel is a 
crucial ingredient in calculating the response functions. A "poorly" behaved spin interaction 
introduces a significantly large model dependence in the neutrino cross sections. 

F , on the other hand, is well behaved and is qualitatively similar to the results of 
microscopic calculations. This is expected, since interactions in this channel are fit to the 
empirical value of the nuclear compressibility K at no- The models SLy4 and SLy5 [Kq — 230 
MeV) were further constrained to reproduce the pressure versus density curve for pure 
neutron matter of the microscopic calculations of Wiringa et al. |27[ . To date, spin dependent 



interactions in the Skyrme models have not been directly constrained. This accounts for 
the large differences between the various models. Whether or not a simple Skyrme model 
can mimic the bulk properties of a microscopic calculation of both spin symmetric and spin 
asymmetric systems remains an open question. For this reason, the use of a Skyrme model 
to describe the spin response of neutron matter for densities of interest in neutron stars does 
not appear to be promising until the spin dependent interaction is constrained. Nevertheless, 
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the energy density functional method provides a consistent framework to study the response 
of hot and dense asymmetric nuclear matter. We hope that further studies will provide 
a better Skyrme energy density functional or an alternative parameterization that will be 
useful for neutron star calculations. 



Microscopic Potential Models 



In a microscopic approach, one starts with the bare interaction and obtains an effective 
interaction by solving iteratively the Bethe-Goldstone equation. Commonly known as the 
Brueckner G-matrix method, this method provides a useful means of arriving at density 
dependent effective p-h interactions. Further refinements include the use of correlated basis 
states and inclusion of a larger class of diagrams In this article, we do not wish to 

address the merits and demerits of these calculations; instead, we contrast the neutron 
matter results obtained by Backmann et al. [^UJ with those of Jackson et al. [23| for the 
Reid-^6 and Bethe-Johnson-f6 potentials (see Fig. ^). In sharp contrast to the Skyrme 
models, these calculations predict a large and repulsive spin interaction (G? > 1). Although 
the qualitative behavior of Fq and M* are similar, important quantitative differences exist. 
In view of this, a knowledge of at least the I = Fermi-liquid parameters of more modern 



calculations pO|,|3Tf would be very useful. 



2. Neutral Current Neutrino Scattering Cross Sections 

In the density regions where spin stability is assured, the neutrino mean free paths may 
be computed for the Skyrme models in degenerate neutron matter, since M*, V 00 , and V w 
are easily expressible in terms of the Skyrme parameters (see Appendix A). The models 
SLy4 and SLy5 have been constrained to better fit systems with large isospin asymmetries. 
For this reason, we choose the model SLy4 to study the neutron matter response. Fig. |] 
shows the neutrino differential cross section for a neutrino energy E u = 25 MeV, at fixed 
momentum transfer q = 25 MeV, as a function of qo/q. At % = 0.04 fm -3 , the p-h 
interaction is attractive in the density-density channel (Fq = —0.7) and repulsive in the 
spin-density channel (Go — 0.6). This accounts for the enhancement (over the Hartree 
result) in the vector channel and suppression in the axial vector channel at small go- The 
enhancement in the spin channel when q /q ~ vp, where vp = kp/M* is the velocity at the 
Fermi surface, is due to the presence of a collective spin excitation that arises due to the 
large and repulsive spin interaction. At finite temperatures, the response is dominated by 
the kinematical region in which go < T. 

The bottom right panels of Fig. ^ show that a strong and attractive interaction in the 
vector channel gives rise to a large peak in the vector response which dominates the total 
differential cross section. The weak collective spin state seen at T=0 is almost completely 
damped at finite temperatures and only modestly enhances the cross sections in this kine- 
matical region. 

The top panels of Fig. |] show results for ng = 0.16 fm -3 . For densities close to no, 
the p-h interaction in both channels is relatively weak. It remains repulsive in the spin 
channel (Gq = 0.15) and attractive in the spin independent channel (F = —0.28). The 
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RPA correlations do not significantly alter the differential cross sections since the modest 
suppression in the spin response is compensated by an enhancement in the density response. 
At higher densities, the Skyrme models investigated here predict a phase transition to a 
ferromagnetic state and become acausal. In our view, these are clear indicators of an inherent 
high density problem in the Skyrme models due to the very strong momentum and density 
dependence of the bare interaction. The large differences in the density dependence of the 
Fermi-liquid parameters for the different models suggest that this is indeed the case. 

In Fig. [5], the neutrino differential cross sections for the neutron matter EOS of Backmann 
et al. [^] are shown. At all relevant densities, the repulsive spin interaction accounts for 
the suppression in the axial response and also for the enhancement in the region where 
qo/q ~ Vp due to the presence of a collective spin resonance. The vector part is enhanced 
at low density and suppressed at higher densities as F Q becomes positive. At 2no, the 
repulsion in this channel is sufficient to produce a weakly damped zero-sound mode. The 
total differential cross sections are significantly suppressed in all cases, and the suppression 
increases with increasing density. 

The collision mean free paths for the Skyrme model SLy4 and for the microscopic model 
of Backmann et al. |29| are compared in Fig. ^| The mean free paths with and without RPA 
corrections are labelled Xrpa and Am*, respectively. The qualitative trends are similar to 
those observed at the level of the differential cross sections. The ratio Xrpa/^m* shown in the 
right panels indicate that typical RPA corrections are large, but model dependent. Neutrino 
mean free paths are suppressed in the model SLy4 due to an attractive p-h interaction. When 
the Fermi-liquid parameters obtained by Backmann et al. |29| are employed, the large and 
repulsive spin interactions enhance the mean free path by as much as a factor of four. The 
differences between the Hartree results (Am*) reflect the different density dependencies of 
the nucleon effective mass predicted by these models (see Fig. ^| and Fig. The density 
dependence of M* in the Skyrme model SLy4 is very strong, but is relatively moderate in the 
microscopic calculation due to Backmann et al. p9] . The increase in the Hartree mean free 
paths with increasing density is a result of the rapidly decreasing nucleon effective mass and 
the use of a nonrelativistic dispersion relation. At high density, a nonrelativistic description 
for the nucleon kinematics breaks down and results in a spurious increase in the neutrino 
mean free path. In the next section, we study the response of matter in relativistic models 
in which some of these deficiencies are remedied. 



B. Relativistic Models 

In Walecka-type relativistic field-theoretical models, nucleons are described as Dirac 
spinors which interact via the exchange of mesons. In the simplest of these models, the a- 
and tu-meson fields simulate the attractive and repulsive character of the nucleon-nucleon 
interaction, respectively. The strengths of the nucleon-meson couplings are fixed by the 
ground state properties of nuclear matter at saturation density by solving the field theory at 



a specified level of approximation ]32] . Further, to account for the isovector character of the 
nucleon-nucleon interaction, the p-meson contribution is included and the pNN coupling is 
fixed by the empirical value of the nuclear symmetry energy. In this section, we study the 
response of neutron matter whose ground state properties are calculated in the mean field 
approximation to the aujp model. 
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The differential cross section per unit volume V for z/-neutron scattering is given by |33 



1 d 3 a G\ E 3 



V d 2 n 3 dE 3 32tt 2 



1 - exp ( -— 



[l-f 3 (E 3 )]xIm(L a ^ ), (10) 



where E% and E 3 are the incoming and outgoing neutrino energies and = Ei — E 3 is 
the energy transfer in the reaction. In terms of the incoming neutrino four momentum 
k = (Ei, k) and the four momentum transfer = (q Q , q), the lepton tensor L a p is given by 

L af3 = 8[2k a k + (k ■ q)g af3 - (k a q p + q a k p ) =p ie^Wq"] . (11) 

The target particle retarded polarization tensor, which is a function of the neutron chemical 
potential, temperature, the kinematical variables qo, and \q\, is 



Imn^ = tanh(J^Jlmn a/3 , (12) 
where Tl a p is the time ordered or causal polarization and is given by 

n Q/3 = -ij y^Tr [T(G 2 (p)J a G i (p + q)Jp)\ . (13) 



The Green's functions Gi(p) (the index i labels particle species) describe the propagation 
of baryons at finite density and temperature. The current operator J M is 7 M for the vector 
current and 7^75 for the axial current. Given the structure of the particle currents, we have 

U aP = c 2 v Ii v ap + c 2 A U% - 2c v c A Il v J . (14) 

For the vector polarization, {J Q , Jp} :: {70,7/3}, for the axial polarization, {J a ,Jp} :: 
{7a 75; 7/375}; an d f° r the mixed part, { J a , Jp} :: {7^75,7/3}. Using vector current con- 
servation and translational invariance, IT^g may be written in terms of two independent 
components. In a frame where = (qo, \q\, 0, 0), we have 



a 2 

n T = Ii v 22 and n L = -f^Ii v 00 

\Q\ 



The axial current correlation function can be written as a vector piece plus a correction 
term: 

K = K + 9^ A - (15) 

The mixed axial vector current correlation function is 

nl u A = ie m0 q a Il VA . (16) 

The above mean field or Hartree polarizations, which characterize the medium response to 
the neutrino, have been explicitly evaluated in previous works j3~i| , |3"o] and are collected in 
Appendix B. At the Hartree level, the only modification that arises is due to the density 
dependent nucleon effective mass. 
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In the aup model, p-h excitations propagate in the medium via the interactions arising 
due to these mesons. The RPA response includes a subclass of these excitations; namely, 
the ring diagrams are summed to all orders by solving the Dyson equation ||. To calculate 
the required polarizations, we begin with the Lagrangian density 

L = Y^Bi-i^d,, - g^Bl^u^ - gpB^K ■ t - M B + g ryB o)B 



+ J2l(-i^d„- mi )l. 



Here, B are the Dirac spinors for baryons and t is the isospin operator. The field strength 
tensors for the ui- and p-mesons are = d^u u — d u uj^ and = d^b u — d u b^, respectively, 
and U(a) represents the scalar self-interactions and is taken to be of the form U(a) = 
(b/3)M n (g aN a) 3 + (c/ '4){g aN aY . The model may then be solved either in the mean field 
approximation (MFT) or in the relativistic Hartree approximation (RHA). In either case, 
the nucleon propagator G{p^) retains the same structure as in the noninteracting case, but 
with constant density dependent shifts to the mass and energy due the presence of scalar 
and vector mean fields. Explicitly, M — > M* = M — g a a and p — > p$ = p — g^uj + \g p b, 
where a, a;, and b are the zeroth components of the a— ,u— , and p-meson mean fields. The 
Hartree polarization functions (see Appendix B) are only sensitive to mass shifts since the 
momentum independent energy shifts exactly cancel. The Dyson equations for the various 
polarizations are 



(17) 



where Tl fl ' u is the RPA polarization function and D Tn is the interaction due to the exchange 
of scalar and vector mesons. The RPA response in these field-theoretical models have been 
studied previously |33|]. Here, we gather the important steps leading to the evaluation of the 
RPA corrections. The correction to the vector part of the response is given by 



n 



n 



(18) 



In order to account for the mixing between the attractive scalar and the repulsive vector 
interactions, the Dyson equation in the isoscalar longitudinal channel becomes 



II s 






tl M 


n 00 _ 





n M n 00 



+ 



ii s fi M 
n M n 00 



-Xs 
-xv 



n M n 00 



(19) 



The interaction matrix is written in terms of the variables 

(W«v) 2 (W m ^) 2 



Xs 



1 - q 2 Jml 



Xu 



1 -Q 2 J m i 



Xp 



{g P /m p y 



(20) 



and xv = I^Xv/l 2 , where xv = Xu + Xp/4 is the sum of u>- and p-meson contributions to 
the vector interaction. Note that the p-meson contribution to the interaction is significantly 
smaller than that of the a;-meson and provides a very small additional repulsion. The 
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dominant contributions arise due to a- and cu-mesons. The solution to the above matrix 
equation yields 

'xv^lo + xs^li - xvxs(n 2 M - n 00 n 5 " 



t2 



e L = 1 + Xs^s + XvRoo + XvXsi^h ~ n 00 n s ) . (21) 

Only the vector mesons contribute to the transverse parts and the RPA corrections in this 
case take the simple form 

<m T = <m 22 = UtXvUt ; e T = l + xv^t ■ (22) 

'/• 

The axial vector polarization is not significantly modified in the aujp model, since the vector 
mesons modify the axial response only through the vector-axial vector mixing. The correc- 
tions to the axial vector response due to the correlations induced by the vector mesons are 
given by the relations 

ft^ = ITJ + IIvXv^va and Il VA = U VA + UvaXv^t- (23) 

To directly modify the axial response, we require a meson with a pseudo-vector coupling to 
the nucleons. The lightest meson which may be coupled to the nucleons in this way is the 
pion. However, the direct part of the pion contribution in this channel will be suppressed in 
the long wavelength limit, since the derivative 7rNN coupling vanishes in the limit —>■ 0. 
In addition, it is well known from nuclear phenomenology that there exists a large repulsive 
component in the spin-isospin channel. To account for this large repulsion, various authors 
(see, for example, Ref. f2~If ) have suggested the use of a very short range repulsive interaction 
parametrized through the so-called Migdal parameter g' . This is usually accomplished by 
modifying the pion propagator according to |36| 



-ir^--r^T-s%- ( 24 ) 

The value of g' is expected to be in the range (0.5-0.7) in nonrelativistic models. However, 
a detailed study to empirically determine this value in relativistic models does not exist. 
For this reason, we introduce a point-like interaction with a pseudo-vector coupling to the 
nucleons and parameterize its strength through the Migdal parameter g'. The strength of 
the coupling is given by xa = g'f^/ m l- The pion contribution, which is expected to be 
attractive albeit small, is not explicitly taken into account, since we are treating g' as a 
parameter that characterizes the total strength in the pseudo-vector channel. With this 
choice, the RPA corrections to the axial vector polarization function is given by 

SU^ = n A T D%Il' A » with D A = (l-D A )- 1 D A and D A = X a9^ ■ (25) 

Having determined the correlated polarization functions, we may calculate the differential 
cross section in a relativistic model. For this purpose, we have chosen a mean field model 
labelled GM3 and a model which includes the vacuum contributions incorporated through 
the relativistic Hartree approximation (RHA). The various coupling constants for the model 
GM3 and a detailed discussion of these field theoretical models can be found in Ref. IfD. 
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1. Neutrino Scattering Cross Sections 



We begin by discussing the results of neutrino mean free paths obtained when RPA 
correlations are ignored (Fig. [j]). In this case, the only medium dependence arises through 
the density dependent M*. With increasing density, the decreasing M* nearly compensates 
for the increase in particle density so as to make the neutrino mean free path approach 
a constant value at high density. This feature, seen in both field-theoretical models, is 
not seen in the free gas case, in which the neutrino mean free paths decrease far more 

2 

rapidly. The nonrelativistic cross sections, which are roughly proportional to M* , are 
reduced significantly with decreasing M*. In the relativistic case, the cross sections are 
proportional to Ep = k\ + M* , resulting in a relatively weaker density dependence. Thus, 
at high densities of interest in neutron stars, the use of nonrelativistic reaction kinematics 
underestimates the cross sections by approximately a factor (M* /Ep) due to the decreasing 
behavior of M* with density. This factor could be as small as 0.5 at high density. 

RPA corrections in models GM3 and RHA due to the a, u, p, and g' correlations suppress 
the cross sections at high density owing to the repulsive character of these interactions. At 
lower density, when M* 3> kp, the mixing between the scalar and vector mesons is large 
and the interaction in the longitudinal vector channel is attractive due to the dominant 
contribution of the scalar meson. At higher densities, mixing becomes unimportant and 
vector meson contributions dominate leading to significant suppression. The transverse 
components of the response in the vector channel are reduced at all densities, since only the 
vector mesons contribute to this correlation function. In Fig. |8| and Fig. [|, the differential 
cross sections with and without correlations are shown for the RHA and GM3 models. In 
the absence of pseudo-vector correlations arising due to g', only the vector part is modified. 

The correction to the axial vector response functions are negligibly small when g' = 0. 
The corrections to the vector part are significant at high density (more than a factor of two). 
However, since the axial current contributions dominate the differential cross sections, even 
large suppressions to the vector part of the response will translate to only modest changes in 
the neutrino scattering rates. To explore the sensitivity to the magnitude of g', we study the 
response for g' = 0.3 and g' = 0.6. Effects due to repulsive correlations are quite significant 
even for g' = 0.3 (see Fig. || and Fig. |9|). 

Fig. [TO shows the mean free paths for the models investigated in this work. The left 
panels show the Hartree cross sections and the right panels show the extent to which RPA 
correlations modify the Hartree results as a function of both density and temperature. RPA 
correlations due to a— ,u— , and p-mesons result in (20 — 30)% increase in the neutrino 
mean free paths over the Hartree results at high density. These modifications are modest 
and show little density dependence, since they predominantly act only on the vector part 
of the neutron matter response. The pseudo-vector corrections modify the dominant axial 
vector response and significantly increase the mean free path. With increasing density, these 
corrections first increase and then decrease due to the small nucleon effective mass at high 
density. 
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III. NEUTRAL CURRENT NEUTRINO CROSS SECTIONS IN 
MULTI-COMPONENT MATTER 



The response of multi-component matter differs in several respects from that of a sin- 
gle component system. At the free gas level, the neutrino cross sections are influenced by 
compositional effects arising due to species specific neutrino-matter couplings, and more 
importantly due to a reduction in the Pauli-blocking. At fixed baryon density and tem- 
perature, additional baryonic components will significantly decrease the degree of particle 
degeneracy and enhance the scattering cross sections. In charge neutral matter containing 
nucleons and leptons in {3— equilibrium, the proton fraction depends very sensitively on the 
dense matter EOS through the density dependent nuclear symmetry energy. The role of 
the nuclear symmetry energy in determining the neutrino scattering and absorption mean 
free paths calculated in the Hartree approximation was investigated in Ref. 0. Here, we 
investigate the effects arising of RPA correlations in a multi-component system. As in Sec. 
II, we study the response in both nonrelativistic and relativistic modes of dense matter. 



A. Nonrelativistic Models 



We begin by considering a two-component system comprised of neutrons and protons. At 
first, we ignore the electromagnetic correlations between the protons and develop the basic 
formalism required to calculate the correlated response of the strongly interacting system to 
neutrino scattering. For a two-component system, the Hartree polarization function may be 
written as a 2 x 2 matrix, diagonal ( U®, n°), where and n° are the proton and neutron 
polarizations, whose real and imaginary parts may be explicitly evaluated using Eq. ([3D and 
Eq. (^), respectively. To compute the 2x2 RPA polarization function, we parameterize the 
p-h interactions in the spin independent and spin dependent channels. The spin independent 
interaction is given by f nn , f pp , and f np , and the spin dependent interaction by g n n,g PP , and 
g np . It must be borne in mind that these parameters are associated with in-medium p-h 
interaction and not the particle-particle interactions. For example, the parameter f nn (g nn ) 
is a measure of the matrix element between two nn~ l states with identical (opposite) spin 
projections, and f np (g np ) measures the matrix element between a nn~ l and pp -1 states 
with identical (opposite) spin projections JJ7 . 



In matrix form, the interactions in the spin 



independent and spin dependent channels are given by 

fpp fpn 



D 



v 



and D , 



gpp gpn 

gpn gnn 



fpn fm 

respectively. The Dyson equation for the vector polarization is 

U + U RPA D n 0_ 



TfRPA 



(26) 



(27) 



The axial polarization is obtained by replacing Dy in the above equation by Da- Solving 
the matrix equations, we obtain 

;i-/ nn n°)n° 

- f n° ) n° 



ttRPA 



1 



f n°n° 

J rip J-J-p 11 n 



f n°n° 

Jpn J-J-p J-J-n 



A, 



f n° - f n° + f n° f n° - f 2 n° n° 

Jnn LL n Jpp p ' J PP p J nn LL n J rip 1L n p 



(28) 
(29) 
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for the vector part, and 



-nRPA 



A. 



1 



[l-g nn U° n )U" p g np U" p Il n 

g pn u p u° n (i- fe n°)n° 



i 9nn n n (^pp rip + g pp rip ^ nn n n n n 



(30) 
(31) 



for the axial part. Multiplying the above polarization matrices by the appropriate vector 
and axial vector couplings for the neutrons and protons, we arrive at the vector and axial 
vector response functions: 



Sv{qo,q) 



YfRPA 



SA{qo,q) 



ttRPA 
11 A 



1 T^Im n£ PA 

- exp {-qo/T) 

(c^) 2 (i - f nn K)K + (c^) 2 (i - f PP n° P )K + 2^f np u° n u 



H' 



1 



— Im TlY A 

- exp {-qo/T) 

(c^) 2 (l - g nn K)K + {c n A )\l - g PP U° p )U° n + 2cf A c n A g np U n U 



(32) 



respectively. In terms of these response functions, the differential cross sections per unit 
volume is given by 



El { l- Um 



(33) 



V dVt 2 dq 8tt 3 

x [(1 + cos9)S v (q ,q) + (3 - cos9)S A (q , q)] , 

where the various kinematical labels appearing above are as in Eq. (fl|). Using Eq. (|33"D, 
the neutrino cross sections may be computed if the p-h interaction is specified. From the 
discussion of the previous section, it is clear that these interactions are model dependent 
and that uncertainties are large. Therefore, we explore different dense matter models to 
identify the generic trends. 



1. Particle-Hole Interactions in Nuclear Matter 



In contrast to neutron-rich matter, for which we have little empirical information, the 
p-h interaction in nuclear matter may be related to empirical values of the Fermi-liquid 
parameters at nuclear saturation density. The p-h interaction is written in terms of the 
Fermi-liquid parameters in Eq. (g). In isospin-symmetric matter, f nn = f pp and g nn = 



This allows us to directly relate the interaction strengths to Fermi-liquid parameters [57]: 



fn 



Fo + F> 



f, 



np 



Fo-F> 



Gp + G'q 



9np 



Gq — G'q 



(34) 



where N = 2M*kf/ir 2 is the density of states at the Fermi surface. Note that we have 
retained only the the I = (s-wave) terms, since we expect the momentum transfer in 
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the p-h channel to be small for low energy neutrino scattering. However, the momentum 
dependent part of the interaction, which gives rise to the I = 1 term, F\ = 3(1 — M*/M) is 
still important since it determines the nucleon effective mass and hence N Q . 

From FLT, we know that the I = Fermi-liquid parameters are directly related to 
macroscopic observables. Fq is related to the isoscalar incompressibility K = 9dP/dnp, and 
Fq is related to the nuclear symmetry energy = (np,/'2)d 2 €/dn\ , where n% = n n — n p is 
the isospin density. Explicitly, 

F = (K/6E F ) - 1 and F$ = (3a 4 /E F ) - 1 , (35) 

with Ep = (k F /2M*) being the Fermi energy. Note that in extracting the p-h interaction 
potential from experimental observables such as K and a^, a consistent value of M* must 
be employed, since both F and Fq depend on M*. 

At nuclear saturation density, investigation of the monopole resonances in nuclei suggest 
that the isoscalar compressibility K = 240 ± 40 Me V ]38| . Information from neutron-rich 
nuclei and observed isovector giant dipole resonances in nuclei require that = 32 ± 5 MeV, 
and empirical determinations of the nucleon effective mass from level density measurements 
in nuclei favor M*/M = 0.7 ± 0.1. For typical values of K = 240 MeV, a 4 = 30 MeV, and 
M*/M = 0.75, the Fermi-liquid parameters are: F = -0.18, Fq = 0.83, and F 1 = 0.75. The 
spin dependent parameter Go, which is related to the experimentally observed isoscalar spin- 
flip resonances is estimated to be small, Go = 0.1 ± 0.1 |39|j40]| . In contrast, the parameter 
G' , which is related to the isovector spin-flip (Giant Gamow- Teller) resonances in nuclei is 
empirically estimated to be large, G' = 1.5 ± 0.2 pT|j39|] . 

Using the empirical values for the Fermi-liquid parameters: F = —0.18, Fq = 0.83, Go = 
0, G' = 1.7, and M*/M = 0.75, we find that the numerical values for the p-h interaction 
strengths in the various channels are: 

fun = fpp = 1-7 x 10" 5 MeV" 2 , f np = -2.7 x 10" 5 MeV" 2 , 

9nn = 9 PP = 4.5 x 10" 5 MeV" 2 , g np = -4.5 x 10~ 5 MeV" 2 . (36) 

Using these values, the neutrino scattering differential cross sections in nuclear matter are 



shown in Fig. [11]. The upper panels are for matter at zero temperature, where only the 
positive go response exists. The suppression due to M* effects are important amounting 
to ~ 50% reduction. Correlations, which are predominantly repulsive in nature due to the 
large empirical values for G' and -Fq, result in significant further reductions. The RPA result 
indicates a factor of two reduction in the differential cross sections in the region q /q vp , 
where v F is the Fermi velocity. The presence of a collective state in the region qo/q ~ v F 
enhances the cross sections in this region. This enhancement, however, is not significant 
enough to override the large suppression seen in the region where qo/q is small. 

Integrating over the q — q space, we obtain the total cross section per unit volume 



or equivalently the inverse collision mean free path. This is shown in Fig. I2|. The left 
panels show the cross sections calculated by taking into account only effects due to M*. 
The results shown are for different temperatures and for a neutrino energy E v = ttT. Right 
panels show the ratio Xrpa/^m*- The resulting increase in A due to the presence of a 
repulsive p-h interaction is approximately a factor of 2.5 at low temperature and decreases 
with increasing temperature. 
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Density Dependence 



The density dependence of the Fermi-liquid parameters is poorly constrained by data. 
Although numerous theoretical models have been constructed to gain insight into their high 
density behavior, there appears to be no general consensus at the present time. Micro- 
scopic calculations of neutron matter differ quantitatively depending on their underlying 
assumptions. These model dependencies are so large that no generic qualitative trends may 
be identified. The exception is the isoscalar parameter F Q , which becomes positive and in- 
creases with increasing density, a feature which may be expected on general grounds as the 
repulsive vector meson contributions dominate. The uncertainties associated with Fq are 
related to the model dependence of the nuclear symmetry energy. In models that favor a less 
than linear increase of 0,4 with density, Fq is expected to decrease with increasing density 
(see Eq. (0)). The state of the art microscopic many-body calculations favor a modest 
increase in the nuclear symmetry energy at intermediate densities ]2"7,3"T]; thus, we may 



expect that Fq will generally decrease. The parameter G' is related to pion condensation, 
since it is a measure the spin-isospin susceptibility of nuclear matter. The large repulsive 
character of G f strongly inhibits s-wave pion condensation in the vicinity of the nuclear 
saturation density. However, at higher densities pion condensation cannot be ruled out a 
priori [|3D|JT|. Thus, while we may expect G' Q to decrease somewhat with increasing density, 



quantitatively it remains very sensitive to the underlying model. The density dependence 
of the isoscalar spin parameter Go, which is not well constrained even at nuclear density, is 
largely unknown. 

Faced with these uncertainties, we begin by assuming that the spin dependent parameters 
are fixed at their empirical values (determined at saturation density), and use schematic 
models to explore the influence of the density dependence of F and Fq. For this purpose, 
we employ a simple parametric form for the EOS (see Appendix A). This model does 
not explicitly address the role of spin dependent interactions and assumes that the favored 
ground state is spin symmetric. In particular, we choose the Skyrme-like models labelled 
'SLn2' with a linear increase in the nuclear symmetry energy. The index 'n' in 'SLn2' takes 
on the values n = 1, 2 and 3 for which K = 120, 180 and 240 MeV, respectively. The 
magnitudes of the RPA corrections to the neutrino mean free paths for these different EOS 
models are shown in Fig. O. Since the dominant contribution to the scattering cross section 



arises from the axial vector response function, the magnitude of the RPA corrections are 
mostly sensitive to the spin dependent parameters. Thus, although the vector response of 
the nuclear medium is modified by about 50-80 % at high density due to RPA effects, the 
changes due to the varying stiffness of the dense matter EOS are small. This suggests that 
the neutrino mean free paths will not be significantly altered due to variations in the nuclear 
compressibility (Fq) or due to variations in the nuclear symmetry energy (Fq) as long as the 
axial contributions are not drastically reduced. 

Fig. [14] shows the behavior of the neutrino mean free paths in symmetric nuclear mat- 
ter for the EOS labelled SL22 as a function of density for the temperatures of interest. A 
comparison of the upper and lower panels show that there is significant temperature depen- 
dence in the Hartree and RPA results, and that RPA corrections decrease with increasing 
temperature. 

To explore the sensitivity to density dependent spin interactions, we employ the Skyrme 
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models discussed earlier. In Fig. [15] the Fermi-liquid parameters in symmetric matter are 
shown. As in the case of pure neutron matter, the results show a large model dependence. 
The top panels in Fig. [16| contain the neutrino cross sections for the model SGII, which has 
been constrained to fit the isovector giant dipole and Gamow- Teller resonances observed 
in nuclei. The neutrino cross sections are enhanced at high density due to an attractive 
interaction in the spin channels. Similar trends may also be expected for the model SkM*, 
since the Fermi-liquid parameters show similar qualitative trends. The models SLy4 and 
SLy5 show distinctly different behavior, especially in the Go and G' Q channels, in addition to 
favoring a larger symmetry energy. Neutrino mean free paths for the model SLy4 are shown 



in the lower panels of Fig. 16. This model is poorly behaved in the spin-isospin channel, 



the negative value for G ! at nuclear density testifying to this fact. Unlike the model SGII, 
which is stable up to ~ 3no (for higher densities Go < — 1), the model SLy4 is stable in the 
spin channel for all densities. However, it becomes unstable to long wavelength spin-isospin 
oscillations (pion condensation) for n > 2n . It is this large attraction that is responsible 
for the significant reductions in the mean free path observed (see lower right panel Fig. |16|) . 

We wish to emphasize that the above model predictions are strictly valid only for those 
densities at which the stability conditions, which require all the i =0 Fermi-liquid parameters 
be greater —1, are satisfied. In view of this, the range of applicability of these models is 
limited to a small window in density centered around the saturation density (see Fig. 15). 
We have investigated a large class of Skyrme models and find this to be a generic feature. 
It is not clear at the present time if the choice of a better parameter set will help alleviate 
this feature. More likely, additional momentum and density dependent terms are required. 
Despite these shortcomings, the Skyrme energy density functional method provides a very 
simple and transparent scheme to calculate both the EOS and consistent opacities required 
for astrophysical problems. We hope that, in the future, a Skyrme parameterization that 
mimics the results of microscopic many-body calculations will become available so as to 
facilitate finite temperature calculations of dense matter at arbitrary asymmetry. 



2. (3— Stable Asymmetric Matter 

We now study the response of charge neutral matter in (3— equilibrium using the 
schematic model labelled SL22, and the Skyrme models SGII and SLy4. The composi- 
tion of matter with and without neutrino trapping, and at finite temperature is calculated 
as outlined in Ref . [I] . The linear increase in the nuclear symmetry energy in the model SL22 
favors relatively large proton fractions compared to those predicted by the SGII and SLy4. 
Using the neutron and proton chemical potentials obtained by solving the EOS, neutrino 
mean free paths are calculated as described earlier. The electrons, which are relativistic 
for all densities of interest, cannot be consistently accounted for in a nonrelativistic treat- 
ment. Therefore, we defer the discussion of the response which includes electrons to the 
next section. 

At first, we ignore effects due to electromagnetic correlations between the charged com- 
ponents and consider only the response of the strongly interacting baryonic components. 
The differential cross sections are computed using Eq. ( P5[ ) with the appropriate baryon 
chemical potentials and effective masses. The p-h interaction potentials in the spin inde- 
pendent channel are computed consistently with the schematic EOS as described earlier, i.e. 
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by taking functional derivatives of the energy density with respect to neutron and proton 
number densities. In general, for asymmetric matter f nn ^ f pp . Therefore, we require three 
independent quantities to characterize the p-h interaction in the spin independent channel 
(unlike only two that were required for nuclear matter) . The three relevant quantities may be 
computed directly, if the form of the potential part of the single particle spectrum is known. 
For the model SL22 the single particle spectrum obtained by a functional differentiation of 
the potential energy density is given in Eq. ( pB| ) of Appendix A. The functional derivatives 
of the single particle potential energy are then related to the p-h interaction parameters as 
follows: 

SUn SU, 5U IL _5U IL 

Jnn- g , hv- 5 » ana Jnp- § - g ■ {<*<) 

For the spin dependent part, we use empirical values since our schematic EOS does not 
explicitly contain spin dependent interactions. These empirical values are determined at 
saturation density and for symmetric nuclear matter. There is a priori no reason to expect 
that they will remain unchanged with increasing density and asymmetry. Our choice is 
therefore primarily motivated by its simplicity, and to a lesser extent by microscopic cal- 
culations that indicate that these parameters do not change appreciably with density |]28 



Fig. [L7| shows the differential cross sections in stellar matter with zero neutrino chemical 
potential. Results for % = n (bottom panels) and for % = 3n (top panels), with pro- 
ton fractions Y p = 0.049 and Y p = 0.22, respectively, are shown. The structure seen in 
the Hartree results (dashed curves) at T = is associated with the fact that the single 
pair component of the the proton response is available only in the kinematical region where 
Qo < \q\ v f, where vf = k F /M* is the velocity at the Fermi surface. For ub = n , the RPA 
cross section is suppressed by roughly a factor of two. The choice of a large and repulsive 
spin interaction is somewhat compensated by the attractive isoscalar component of the p-h 
interaction. For % = 3n , p-h interactions are repulsive in all channels and account for the 
larger (factor of four) suppression seen in the upper panels of Fig. [17, 

In the Skyrme models, the p-h interaction in all the relevant channels is explicitly cal- 
culable as functions of density and proton fraction, and are given in Appendix A. The 
differential cross sections are computed using Eq. ( |33| ) with the appropriate model depen- 
dent inputs, i.e. the p-h interaction, the neutron and proton chemical potentials and the 
effective masses. Results for the models SGII and SLy4 are shown in Fig [IB]. As a result 
of the large density and asymmetry dependence of the p-h interaction predicted by Skyrme 
models, the RPA corrections vary significantly with changing density and proton fraction. 
The top panels show results for tie = 0.16 fm -3 and proton fraction Y p = 0.05. The dif- 
ferences seen in the Hartree results are due to the different effective masses and proton 
fractions predicted by these models. The RPA results differ qualitatively between these two 
models. At nuclear saturation density (left panels), both models are stable with respect 
to spin-isospin fluctuations and contain nearly equal proton fractions. The RPA response 
differs mainly due to the different spin-dependent interactions: SLy4 favors a large and 
attractive spin-isospin interaction, which strongly correlates the neutron and proton spins 
and accounts for the enhancement in the kinematical region where the protons contribute 
(region where qo/q < v F ). In model SGII, the situation is reversed since the p-h interaction 
is weakly repulsive in the spin channels and accounts for the modest suppression seen in the 



left-lower panel of Fig. 18. At 2hq, the p-h interaction is large and attractive in both the 
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isoscalar and isovector spin channels for both models. For model SLy4, this turns out to be 
so large that the spin-symmetric j3— equilibrium ground state is unstable to long wavelength 
spin-isospin fluctuations, a feature which was already encountered in the case of symmetric 
matter. Although the spin forces are large and attractive even for the model SGII, they are 
not sufficiently large to energetically favor a spin-asymmetric state. The resulting enhance- 
ment in the RPA response is about a factor of 4 at small go (note that for T=10 MeV, the 
RPA response plotted has been multiplied by 0.25) and decreases with increasing go. 

The collision mean free paths are shown in Fig. [19]. The Hartree mean free paths are 
shown in the left panels and the ratio of the RPA mean free paths to the Hartree results 
are shown in the right panels. The results show trends similar to those observed earlier for 
neutron matter and for symmetric nuclear matter. Results for the Skyrme models are not 
shown, since the RPA modes become unstable with increasing density and asymmetry. At 
low density and for small asymmetries, the trends are similar to those shown in Fig. 



B. Relativistic Models 

The basic formalism required to calculate the response of a system containing nucleons 
and electrons within the framework of a relativistic mean field theory may be found in Ref. 
H. Here, we summarize the main results of their work and extend it to incorporate spin 
dependent correlations. 

Neutron and proton single pair excitations are coupled due to strong interactions. In 
addition, the proton and electron p-h excitations are coupled due to electromagnetic inter- 
actions. To incorporate these correlations consistently, we write the various polarizations 
functions in matrix form. The vector part of the polarization has transverse and longitudinal 
components. The transverse part is affected only due to the vector mesons. For a system 
comprised of neutrons, protons, and electrons, the RPA polarization is given by the Dyson 
equation 

u PPA = n T + nf PA d v n T , (38) 

where the three-component transverse polarization matrix is diagona^n^, n^, n^) and the 
interaction matrix, defined in terms of x 7 = — e 2 /g^, xv — Xu> + Xpi an d Xi — Xu> — Xpi is 
given by 

— X7 X7 



X-y -Xy ~ Xv -Xi 
-Xi ~Xv 



(39) 



The above equation accounts for electromagnetic, isoscalar, and isovector strong interaction 
correlations. Since the scalar a- and the vector a;-mesons mix strongly in the longitudinal 
channel, the polarization matrix is written as a 4 x 4 matrix equation |4T|. Explicitly, the 
Dyson equation for the longitudinal part is 

Hf 5 A = H L + Hf PA D L H L , (40) 

with the longitudinal polarization for the n,p, and e system given by 
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(41) 



and the interaction matrix 



D, 



-xs 

X 7 







X 7 - Xv -Xi 
-Xi ~Xv 



(42) 



The interaction parameters are: xv = qftXv/q 2 , Xi — iftXi/l 2 , an d X 7 = I^X-f/l 2 - The 
dominant modification to the axial polarizations at small momentum transfers arises due to 
the pseudo-vector contact interaction discussed earlier and is parameterized by the Migdal 
parameter g' . Further, since the electromagnetic contributions to spin correlations are neg- 
ligible, we account only for the correlations due to strong interactions between neutrons 
and protons. It is conventional to introduce the short range correlations by modifying the 
pion propagator in the pseudo-vector form of the pion-nucleon interaction described by the 
Lagrangian density 



L = -(f*/m n )ilry s 'fTil>d IM 'ir . 
The Migdal parameter g' enters through the modified pion propagator 



q m u z -m^+te 
and the NNtt vertex becomes 



1 NNtt 



q m W 
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V / 2(A/m^) 75 7^ 



(43) 



(44) 



(45) 



At small momentum transfers, the repulsive contact term overwhelms the attraction arising 
due to one-pion exchange. For this reason, we retain only the contribution from the g' term 
in the RPA equations. For a two-component system, this requires us to extend the Dyson 
equation for the axial vector part in Eq. ([25]) to an 8 x 8 matrix. However, it can be explicitly 
shown that for (? A = the cross terms cancel and the RPA axial polarization for the n-p 
system closely resembles that for a one-component system, but with = 11^ + IT^. This 
greatly simplifies the formalism required to incorporate spin dependent correlations. The 
RPA axial polarization is therefore once again given in terms of Eq. ( |25| ) . 

We have ignored effects due to mixing of the vector and axial correlations, since they are 
proportional to Hva, which is negligible at small momentum and energy transfers. 

Substituting the vector and axial vector RPA polarizations for the correlated 3 compo- 
nent system into Eq. (|T0D, the differential cross section for the npe system is 
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1 - exp(-q /T) 



(46) 
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Utilizing the particle fractions, nucleon effective masses, and particle chemical potentials in 
(3— equilibrated stellar matter, the neutrino cross sections and scattering mean free paths 
are straightforwardly computed using the above equation. 

In the model labelled RHA, the effects arising due to vacuum loops are explicitly taken 
into account both at the level of the EOS and the response functions. The differential cross 
sections for this model are shown in Fig. ^0|for the model GM3, wherein nonlinear scalar self- 
interactions play an important role pB| (it acts to lower the compression modulus, increase 
the nucleon effective mass and make the "effective" cx-meson mass density dependent). Since 
the Migdal parameter g' plays a dominant role in determining the RPA response, results 
for different values of g' are shown. The modification to the vector response due to scalar, 
vector, and isovector mesons is large, but, it does not translate to large suppressions in the 
differential cross sections (top left panel of Fig. ^0|), since the response is dominated by the 
axial current coupling of the neutrinos. The total suppression remains very sensitive to g' for 
which empirical values are expected to be in the range 0.5 — 0.7. RPA correlations suppress 
the neutrino scattering cross sections by typically a factor of 2 — 3. 

The RPA and Hartree mean free paths of thermal neutrinos for the model GM3 are 
compared in Fig. [H]. Effects due to correlations increase with density and decrease with 
temperature. The decrease in the ratio Xrpa/ (lower panels) at high density is due to the 
decreasing M*. The pseudovector point-like coupling introduced through g' is sensitive to 
M* [36] and its effects decrease with increasing density. The results in the lower left panels 
indicate that the RPA correlations, even for a large g', will typically increase the mean free 
path by a factor 2 — 2.5. 



IV. CHARGED CURRENT NEUTRINO CROSS SECTIONS IN 
MULTI-COMPONENT MATTER 

The discussion thus far has fo cussed upon the role of correlation effects on the neutral 
current scattering reactions. While these reactions are the only source of opacity for the \x 
and r neutrinos, the electron neutrino mean free path is for the most part dominated by the 
charged current reaction v + n — > e~ + p. The absorption reaction is kinematically different 
from scattering, since the energy and momentum transfers are not limited by the matter's 
temperature alone. The energy transfer is typically of order fx = \x n — \x v . In the extremely 
degenerate case, for which T j \i << 1, and in neutrino-poor matter, final state Pauli-blocking 
due to electrons and momentum conservation restricts the available phase space. During 
the deleptonization epoch, the typical neutrino momenta are large (~ 100 — 200 MeV) and 
despite the mismatch of the neutron, proton, and electron Fermi momenta, momentum 
conservation is easily satisfied with the aid of the neutrino's momentum. This is no longer 
true in the late stages of the cooling phase, in which the neutrino energies are of order ksT. 
Further, a finite neutrino chemical potential affects the relative displacement between the 
neutron and proton chemical potentials according to fx = jx e — jx u . 
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A. Non-Relativistic Models 



We begin by extending the formalism developed in Sec. Ill to incorporate effects due to 
correlations. The differential cross section per unit volume, for the absorption reaction in 
matter containing neutrons, protons, and electrons follows from Eq. (]]]), and is given by 
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g v (1 + cos 9)S v (q , q) + g 2 A (3 - cos 0)S A (q o , q) 



(47) 



where gy = 1 and g A = 1-26 are the charged current vector and axial vector coupling 
constants. The particle labels 2 and 4 correspond to the neutron and proton, respectively, 
and 3 to final state electron. The other kinematical variables are similar to the case of 
scattering. The isospin and spin-isospin density response functions are defined in terms of 
the asymmetric polarization IF. Due to the difference between neutron and proton chemical 
potentials, and of the single particle energies, IT differs significantly from II . For the simple 
form of the single particle energy given by Eq. (Q), an analytic expression for the imaginary 
part may be derived (see Ref. |] for details), 
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(50) 



arise due to the isospin dependence of the single particle potential energies. The real part 
is calculated numerically by using the Kramers-Kronig relation Eq. (|7|). The RPA response 
functions are related to the Hartree polarization IT by the Dyson equation. Explicitly, 
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(51) 
(52) 



where the p-h interaction in the vector and axial vector channels are given by f c and g c , 
respectively. 



Symmetric Nuclear Matter 

To elucidate the role of correlations, we begin with a brief discussion of the neutrino 
absorption reactions in symmetric nuclear matter. Here, the kinematical complications 
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arising due to dissimilar neutron and proton Fermi surfaces and electron blocking are absent. 
To begin, we once again need the strength of the p-h interaction for the spin independent 
and the spin dependent channels. In contrast to scattering, both charge or isospin are 
transferred along the p-h channel in the charged current absorption reaction. Thus, only 
the isospin changing part of the p-h interaction is relevant for the charged current reaction. 
In terms of the Fermi-liquid parameters, this implies that only Fq and G' contribute. For 
zero range interactions, the potentials required to calculate the RPA response are given by 
f c = 2Fq/Nq and g c = 2G' /N , where the factor two arises due to isospin considerations. 
In Fig. the differential cross section for the charged current neutrino absorption reaction 
in symmetric nuclear matter is shown. As discussed earlier, large and positive empirical 
Fermi-liquid parameters indicate a correspondingly large and repulsive p-h interaction. This 
accounts for the large suppression seen at small go- The suppression is roughly proportional 
to (1 + Fq)~ 2 in the vector response and (1 + G' )~ 2 in the axial response. At high q and 
for small temperature, well defined collective modes in the isospin-density and spin-isospin 
density channels appear. These correspond to the well known Giant Dipole (labelled as 
GD in Fig. |22[) and Giant Gamow- Teller (labelled as GT) resonances found in nuclei. With 
increasing temperature, these collective states are significantly broadened due to Landau 
damping (see the right panel of Fig. ^) . 



Asymmetric Nuclear Matter 

In asymmetric nuclear matter, strong interactions affect the charged current opacity both 
because of correlations as discussed above, and because of the important role they play in 
determining the quantity ft = \i n — fi p . The latter, was discussed earlier in Ref. |l|]. In 
the following discussion, we show that the p-h interaction in the vector channel and fi are 
directly related. It was noted in Eq. ( |35|) that Fq was related to the nuclear symmetry 
energy. The nuclear symmetry energy in turn determines fi, and thereby the proton fraction 



for stellar matter in (3— equilibrium [5T]. Thus, a particular choice of the p-h interaction in 
the isospin-density channel uniquely determines both fi and the extent to which correlations 
suppress the cross sections. Relative to the case in which no interactions are considered, these 
effects oppose each other. From Eq. (^), it may be verified that fi is linearly proportional 
to (1 + Fq). Correlations suppress the isospin density fluctuations by a factor roughly 
proportional to (1 + Fq)~ 2 . The cumulative effect on the total absorption cross section is 
easily deduced by noting that the cross section is a linear function of fi (see Ref. |]J), and is 
inversely proportional to (1 + Fq) 2 . This underscores the need to employ a p-h interaction 
that is consistent with the EOS employed to compute the composition of charge neutral 
(3— equilibrated matter. Here, we employ the schematic potential model SL22 to compute 
the composition and the p-h interaction in the spin independent channel. 

The p-h interaction in the spin-isospin channel has been studied extensively in the con- 
text of charge-exchange nuclear reactions |39|, muon capture rates on nuclei p2| , and pion 



condensation In these studies, a p-h interaction arising due to the exchange of tt- and 
p-mesons are included to account for the momentum dependence, in addition to the short 
range repulsion parametrized through the Migdal parameter g'. The isovector interaction 
in the longitudinal channel arises due to it exchange and in the transverse channel due to 
p-meson exchange. We employ this widely-used form for the p-h interaction, which has been 
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successful in describing a variety of nuclear phenomena |44j] . The longitudinal and transverse 
potentials in the it + p + g' model are given by 

v ^=§X^h*/ M+e ) (53) 

VTM = ^l 4-7-™f > {q) + 9 )> (54) 

where C p = 2, F„ = (A 2 - m 2 )/(A 2 - g 2 ) and F p = (A 2 - m 2 )/(A 2 - g 2 ) are the ttNN and 
pNN form factors |44[], and the numerical value of g' used was 0.6. With this choice for the 
p-h interaction, the axial response function is given by 



Sa(Qo, q) 
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ImIF(g ,g)(^- + J- ) (55) 



1 -exp[-(g + /i)/T] 
e L = [1 - 2K L Re IT (g o , g)] 2 + [2K L Im IT (g o , q)f , (56) 
e r = [1 - 2V^Re IT (g o , g)] 2 + [2VrIm n (g o , g)] 2 . (57) 

Using the axial response function given above and a p-h interaction in the isospin-density 
channel, which is directly related to the nuclear symmetry energy predicted by the EOS, 
the neutrino absorption cross sections are computed. 

In Fig. ^3], the Hartree and RPA results for the differential cross section for the absorption 
reaction are shown for ris = 0.16 fm -3 and = 0.48 fm -3 . The neutrino chemical potential 
and proton fraction determined by the finite temperature EOS are also shown in the figure. 
Since only those neutrinos close to the Fermi surface participate in transport, we have chosen 
the neutrino energy E u = p u . The response peaks in the region where go ~ ~f L - The RPA 
suppression is roughly a factor of two for the kinematics shown here. The relatively smaller 
suppression seen here is due to the large energy and momentum transferred to the baryons. 
Since the charged current probes small distances, the effects of many-body correlations are 
somewhat suppressed. Further, the pion exchange contribution to the p-h interaction is 
attractive at large q and acts to decrease the short range repulsion due to g' . 

The neutrino absorption mean free path is shown in Fig. ^3] for neutrino trapped matter 
at Yl = 0.4 for different temperatures and densities of relevance. The trends are very similar 
to the case of scattering. The left panel shows the Hartree results, wherein effects arising due 
to the density dependent nucleon effective masses and single particle potentials are included. 
The increase in the mean free path at high density arises due to the rapidly dropping M*(ub) 
and the use of a nonrelativistic form for the single particle energy. The magnitude of the 
RPA corrections are shown in the right panels. On average, the enhancement is about 
a factor 2 — 3. At intermediate baryon density, this enhancement increases with density 
and decreases with temperature. At higher density, the rapid decrease in M* decreases the 
density of p-h states at the Fermi surface and thereby reduces the magnitude of the RPA 
corrections. 

During the cooling epoch of a protoneutron star, the neutrino chemical potential becomes 
progressively smaller. For typical electron neutrino energies of order T, the kinematical re- 
striction on the charged current process becomes important at low temperature. The density 
dependence of the symmetry energy essentially determines the density and temperature at 
which these kinematical restrictions severely inhibit the charged current rates. In the model 
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SL22, which is characterized by a linear symmetry energy, the kinematical restrictions are 
unimportant for % > 2. At lower densities, and for T <C fi p , there is significant suppression. 

In Fig. |25], the Hartree and RPA differential cross sections are compared. The momentum 
transfer was held constant at \q\ — fi e + nT, to ensure momentum conservation. The response 
peaks in the region q ~ (— fi e + E u ) and exponentially decreases for larger q due to Pauli- 
blocking of the degenerate electrons. On the other hand, the final state Pauli-blocking due 
to protons is negligible, since protons are fairly nondegenerate. This, combined with the 
large negative energy transfers, results in the real part IT becoming large and negative, 
leading to significant reductions (factor ~ 4) in the RPA cross sections. For = 0.48 
fm -3 , a larger proton fraction and a lower nucleon effective mass decreases the magnitude 
of the real part IF. In this case, the suppression due to correlations is only about a factor 
of 1.5 - 2. 

In Fig. neutrino absorption mean free paths for matter in which {i v = are shown. 
The trend identified in the preceding discussion, by which we may expect to find a large 
enhancement at intermediate densities and low temperatures, is clearly seen. For n# > 2n Q , 
correlations enhance the mean free paths by a factor ~ 3 — 4. For intermediate densities, 
the enhancement is significantly larger and very sensitive to the ambient temperature. 



B. Relativistic Models 

In field theoretical models, isovector correlations are induced by interactions mediated 
by the p-meson. Axial vector correlations, in the spin-isospin channel, arise due to the 
presence of a strong short ranged component, which is commonly parameterized as a contact 
interaction of strength g' . The calculation of these vector and axial vector responses is 
similar to that discussed in Sec. IIB, but with two important differences. First, the a- 
and o;-mesons do not contribute when isospin is transferred along the p-h channel. Second, 
the polarization functions now describe the propagation of neutron-hole and proton-particle 
states, which depend on fi = /i n — /i p and the isovector part of the single particle potential 
energy. The relevant polarizations are collected in Appendix B. 

The differential cross section for the charged current absorption reaction is the same as 
that for neutral current scattering, Eq. (0), with the modification that f3(E 3 ) is replaced 
by f e (E e ). The vector part of the correlation function Efig is modified due to the p-meson 
contribution. Explicitly, the correction to the longitudinal part is 

= -=g|Sg (58) 
1 + 2x/IIoo 

and to the transverse part is 

ATT att 2II r xjII r 

SUt = 51122 = i + 2 X ,n r • (59) 

For momentum transfers of interest, the axial vector current correlation function is modified 
mainly by the large and repulsive contact interaction of strength g'. As discussed earlier, 
this is conventionally introduced as a correction to the pion propogator (see Eq. (|24])). In 
principle, the contact term will lead to p— tt mixing via in- medium particle- hole states. These 
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mixing effects are ignored in this work. This is justified both because g' is a phenomenological 
parameter, which is determined by ignoring these mixing effects, and because the mixing, 
which is proportional to the vector-axial vector polarization function, is negligible at small 
momentum transfers. The RPA equation for the axial vector polarization is as in Eq. ( p5|) . 
but with xa replaced by 2xa, where the factor 2 arises from isospin considerations. 

In Fig. |27|, the RPA and Hartree differential cross sections for neutrino absorption are 
shown for ris = no and Ub — 3no, respectively. The results are presented for lepton rich, 
neutrino trapped matter with Yl = 0.4. The electron neutrino energy and the momentum 
transfer are both set equal to the neutrino chemical potential. The composition, individual 
chemical potentials, and the density dependent baryon effective masses are computed us- 
ing the field theoretical model GM3. The RPA differential cross sections are significantly 
suppressed relative to the Hartree results, particularly around the peak position i.e., when 
% ~ — (a*ti — /%>)■ The suppression decreases with increasing density. 

The absolute magnitude of the differential cross section depends sensitively on the neu- 
trino energy, the electron chemical potential, composition, and temperature. In order to 
assess the overall effects of RPA correlations, both as a function of density and temperature, 
we integrate over qo and q and compare the RPA and Hartree results for the neutrino mean 
free paths in neutrino-trapped matter in Fig. ESI. As expected, RPA correlations increase 



the neutrino mean free paths. The trends are very similar to those seen in the nonrelativis- 
tic models. This is mainly due to the fact that g' plays a dominant role in suppressing the 
axial vector response in both cases. At low density, the minor differences seen between the 
two approaches may be attributed to the different p-h interactions employed. Relativistic 
effects become important with increasing density and this accounts for the differences seen 
at higher density. 

The RPA and Hartree differential cross sections in neutrino poor matter (fi u = 0) are 
shown in Fig. |29. Due to important compositional and kinematical differences, the neutrino 
absorption cross sections differ from those of nonrelativistic models. At low density, Pauli 
blocking due to protons is small due to their low concentrations, which acts to enhance the 
polarization of the medium in the isospin channel. The suppression at % = n (left panel) 
is approximately a factor of 4. The suppression is reduced at high density; at n B = 3n 
(right panel), the suppression is approximately a factor of 2. 

The RPA and Hartree absorption mean free paths are compared in Fig. |30] as a function of 
density and temperature. The qualitative trends are again very similar to those encountered 
in the nonrelativistic case (see Fig. Ej]). The quantitative differences are due to the different 
compositions arising due to the different isovector interactions employed. It is clear that 
relativistic effects play an increasingly important role with increasing density, decreasing the 
level of suppression compared with the non-relativistic case. 



V. EFFECTS ON PROTONEUTRON STAR COOLING 

Within a tenth of a second of the core bounce which precedes a supernova explosion in the 
death of a massive star, the newly- formed neutron star, or protoneutron star (PNS), becomes 
quasistatic since its structure further evolves only in response to thermal and compositional 
changes in its interior. These changes are, for the most part, controlled by the diffusion of 
neutrinos from the interior. At high density and for temperatures above several MeV, the 
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neutrino-matter reaction rates are sufficiently high to ensure that neutrinos are in thermal 
and chemical equilibrium. Thus, in the equilibrium diffusion aproximation (EDA), the 
neutrino distribution function in these regions is both nearly Fermi-Dirac and isotropic. 



These assumptions have been used to study the evolution of the PNS 

In the EDA, the distribution function in a spherically symmetric background can be 
approximated by an expansion in terms of Legendre polynomials to 0(/i), with 11 being the 
angle between the neutrino momentum and the radial direction. Explicitly, 



/(<*>, ^) = /oM +/i/iH , fo = (1 + exp[(u - n v )/T\) 
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(60) 



where fo is the Fermi-Dirac distribution function at equilibrium (T = T mat , fi u = nl q ), and ui 
and \x v are the neutrino energy and chemical potential, respectively. This allows the energy- 
integrated fluxes of lepton number and energy to be expressed as a linear combination of 
gradients of the degeneracy parameter 77 = fi u /T and the temperature 
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The diffusion coefficients D2, -D3, and -D4 arise naturally from the transport equations and 
contain all the microphysics of the neutrino-matter interactions. They are defined by 



dx ^( w )/oM(i-/oM) 



(63) 



where x = uo/T and D(oS) = {n a + k s )~ 1 . The latter term includes the contribution of both 
absorption-emission opacities (K a ) and scattering opacities (k s ), which are related to the 
relevant cross sections by 



l-/o(a;) V 



(64) 



The fluxes in Eqs. (65) and (66) are technically for one particle species. To include all six 
neutrino types, we redefine them in an obvious notation as 



D 2 = m + , D 3 = - D^C , D A = Dl* + + AD"/ . 



(65) 



Effects on PNS evolution produced by changes in the neutrino opacities can thus be under- 
stood in terms of changes to the diffusion coefficients. 

In Fig. |3l] we show the ratio of each of the diffusion coefficients calculated in the RPA 
approximation to those calculated in the Hartree approximation. The left panels display 
the results for neutrino-free matter and the right panels show the results for neutrino-rich 
matter (Y u = 0.06). In the region of interest during the first, or deleptonization, phase 
of the PNS evolution, in which the trapped neutrino abundances are high and 10 MeV 
< T < 50 MeV, the RPA corrections to the diffusion coefficients are in the range 1.3- 
1.75 for matter above nuclear density. This correction, at first glance, seems smaller than 
expected; a correction factor w 2 was obtained for the mean free paths at E u = ttT (see 
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Fig 21). However, the RPA correction redistributes strength in the spectrum (see Fig. 20), 
decreasing the cross section at energies near the peak while increasing the cross section for 
larger energy exchange. The diffusion coefficients result from an energy integration over the 
entire spectrum. Another important feature is that in the region below about half nuclear 
density, and at low temperatures, the RPA corrections are diminished since the composition 
of matter becomes dominated by nuclei instead of the npe liquid. Therefore, the diffusion 
coefficients in the outer regions of the PNS will not be strongly altered by RPA corrections. 

During the second, or cooling, phase of PNS evolution, in which the excess neutrinos 
have been largely lost, the neutrino chemical potential is relatively low and nearly constant 
throughout the star. The evolution is then governed by the neutrino energy flux, Eq. (|62|). 
in which the electron neutrino chemical potential gradient can be ignored. Thus, D± will 
become the dominant diffusion coefficient. The temperatures are now smaller than in the 
earlier, deleptonization phase, being in the range 1 < T/MeV < 20. Fig. |3l| shows that 
the RPA corrections, at three times nuclear density (0.5 fm -3 ), are larger than in neutrino- 
rich matter, and can exceed a factor of 2.5. Nevertheless, at and below nuclear density, 
the corrections are considerably smaller, as observed in the previous case. The effect of 
RPA corrections may be expected to be more evident in the cooling phase than in the 
deleptonization phase. 

To study the effect in a realistic PNS simulation, we approximated the functions 
Df p / Df , for i = 2 — 4, by a fit in terms of n^, T and Y v . Fig. |32] shows the func- 
tion assumed for Y v = 0; the fits for other values of Y v are similar. Notice the strong density 
dependence in the region ris < n . We used these fits to correct the diffusion coefficients 
calculated in Ref. and performed a new PNS cooling simulation with the same code. 



We also wished to compare these results to those of Burrows and Sawyer |T2[, who also 
studied the effects of correlations in the RPA approximation but simply assumed that the 
total opacities were decreased by a factor of 3.33 above either 1/2.7 or 1/5.4 times nuclear 
density, respectively. Therefore, we performed two additional exploratory simulations with 



the identical assumptions. The results of all three simulations are displayed in Figs. 33 - 35 



In Fig. |33], the evolution of the central values of the important thermodynamic quantities 
is shown. The solid line corresponds to the baseline model, the dashed and dot-dashed lines 
to the exploratory models A and B, and the dotted line to our RPA calculation. The main 
effects of the larger mean free paths produced by the exploratory models and by our RPA 
corrections are that the inner core deleptonizes more quickly. In turn, the maxima in central 
temperature and entropy are reached on shorter timescales. In addition, the faster increase 
in thermal pressure in the core slows the compression associated with the deleptonization 
stage, although after 10 s the net compressions of all models converge. 

The relatively large, early, changes in the central thermodynamic variables do not, how- 
ever, translate into similarly large effects on observables such as the total neutrino luminosity 
and the average radiated neutrino energy, relative to the baseline simulation. The luminosi- 
ties for the different models are shown as a function of time in Fig. |34| and Fig. Figure |4 
shows the early time development in detail. The exploratory models agree with the results 
reported by Burrows and Sawyer |12| . However, the magnitude of the effects when our full 
RPA corrections are applied is somewhat reduced compared to the exploratory models. It 
is especially important that at and below nuclear density, the corrections due to correla- 
tions are relatively small. Since information from the inner core is transmitted only by the 
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neutrinos, the time scale to propagate any high density effect to the neutrinosphere is the 
neutrino diffusion time scale. Since the neutrinosphere is at a density approximately 1/100 
of nuclear density, and large correlation corrections occur only above 1/3 nuclear density 
where nuclei disappear, we find that correlation corrections calculated here have an effect 
at the neutrinosphere only after 1.5 s. Moreover, the RPA suppresion we have calculated 
is considerably smaller than those reported in Ref. [12], reaching a maximum of about 30% 



after 5 s, compared to a luminosity increase of 50% after only 2 s. However, the corrections 
are still very important during the longer-term cooling stage (see Fig. |35|), and result in a 
more rapid onset of neutrino transparency compared to the Hartree results. 



VI. CONCLUSIONS AND OUTLOOK 

The effects due to correlations are important, albeit model dependent. Due to the 
large neutrino-baryon axial vector coupling, the residual interaction in the spin dependent 
channels plays a crucial role in determining the magnitude of the RPA corrections. In pure 
neutron matter, these model dependencies were carefully studied. In the Skyrme model, 
the p-h interaction in the spin channel becomes large and attractive with increasing density. 
This drastically reduces the neutrino mean free path by enhancing the spin response. The 
opposite behavior is suggested by variational calculations of the Fermi-liquid parameters in 
neutron matter. In this case, the residual interaction at high density is strongly repulsive in 
the spin channels. A repulsive p-h interaction generically shifts the strength from regions of 
small energy transfer go to large energy transfer. In the case of strong repulsion, collective 
spin excitations arise and the strength at small go is highly suppressed. 

Not unexpectedly, the Skyrme models, in which the force parameters are calibrated only 
at nuclear density, either become unstable to spin fluctuations or violate causality at high 
density. Schematic models such as SL22, which are partly based on the Skyrme approach, 
are somewhat better behaved at high density, but they do not address the issue of spin 
dependence in the nucleon-nucleon interaction. In these models, the RPA corrections to 
the vector response may be calculated consistently with the EOS. For the axial part, an 
empirical spin-isospin interaction determined by the value of the Fermi-liquid parameter G 
in symmetric nuclear matter may be used. 

The situation is quite similar in field-theoretical models, wherein the full relativistic 
structure of the baryon currents and more importantly, relativistic kinematics of the scat- 
tering reactions are naturally incorporated. In these models, the scalar-attractive contribu- 
tion dominates at low density leading to enhancement in the vector part of the scattering 
response. With increasing density, the scalar decouples and the repulsive vector mesons 
suppress the cross sections. This is a generic feature. The RPA corrections in the vector 
channels generally suppress the neutrino cross sections at high density and moderately en- 
hance it at low density. In addition, the RPA corrections in the vector channel and the EOS 
properties such as the compressibility and nuclear symmetry energy are intimately related. 
While this stresses the need for consistency between the EOS and neutrino opacities, there 
appears to be some degree of decoupling between the strong interaction effects that influ- 
ence the EOS and those that influence the neutrino opacities. This is almost entirely due 
to the different roles played by a repulsive spin dependent interaction in EOS and opac- 
ity calculations. They directly influence the dominant part of the neutrino opacities but 
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play little role in determining the properties of spin symmetric matter as long as the spin 
interaction is repulsive and the spin correlation energy is negligible compared to the other 
interaction energies. This is true for the field-theoretical model, wherein the introduction of 
a strong repulsive pseudovector-isovector interaction through the Migdal parameter does not 
contribute to the EOS at the mean field level, but plays an important role in determining 
the neutrino mean free path. 

The RPA correlations in multi-component matter are qualitatively different from those 
in single-component system. The presence of even a small proton fraction affects the shape 
of the matters response function. The generic trend is that it increases the strength at small 
energy transfer. The inclusion of electromagnetic correlations does not significantly affect 
the neutrino mean free paths. The mixing of the photon with the vector mesons via the p-h 
states could result in the presence of plasmon like collective states in the system. Although 
this might play an important role in determining the average energy transfer in neutrino 
matter interactions, it does not affect the overall neutrino mean free paths. 

The dominant role played by spin dependent nucleon-nucleon interactions in the medium 
and its present poorly constrained status does not allow us to draw definite conclusions. For 
the vector part, we showed that RPA opacities are consistent with the underlying EOS. In 
subsequent works, we hope that better constraints on the spin dependent interactions in 
dense matter will allow calculations which will pin down its effects on both the EOS and on 
the opacities. At the present time, however, the RPA correlations remain sensitive to this 
poorly known density and isospin dependence of the p-h spin interaction. 

Notwithstanding the above caveats, and the differences between nonrelativistic and rel- 
ativists approaches such as the spin- and isospin-dependent interactions and the nucleon 
effective masses, suppressions of order 2-3, relative to the case in which correlations are 
ignored, are obtained. It is also satisfying that the independent investigations of Burrows 
and Sawyer |12| have led to a similar conclusion. 

We studied the effects of RPA correlations on the evolution of protoneutron stars. Our 
calculations indicate that while RPA correlations suppress the neutrino cross section by a 
factor « 3 at thermal neutrino energies, the net effect on energy-integrated quantities such 
as neutrino diffusion coefficients is more moderate. Moreover, the effects of correlations are 
diminished below nuclear densities, and disappear in low-density matter that is dominated 
by nuclei. Thus, we observe significant effects in observable quantities such as neutrino lumi- 
nosities only after 1 second, and it is only after several seconds that the neutrino luminosities 
are increased by as much as 30%. Larger effects occur during the late-time thermal cooling 
phase, and the protoneutron star will become transparent on shorter timescales when RPA 
corrections are included. 
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APPENDIX A: PARTICLE-HOLE INTERACTIONS 



The p-h interaction in the Skyrme and Skyrme-like schematic models may be directly 
obtained by double functional differentiation of the potential energy density. The potential 
energy density for the density and temperatures of interest is mainly a function of baryon 
density and proton fraction. The temperature dependence, which enters only through the 
explicitly momentum dependent interactions, is weak. In what follows, we provide analytical 
expressions for the p-h interaction at zero temperature. 

The effective nucleon-nucleon interaction in the standard Skyrme model |23| is given 
in Eq. (|9|). The potential energy density in the Hartree-Fock approximation for the above 
interaction may be computed using the standard method by employing plane wave states 
for the nucleons [23]. In terms of the neutron and proton Fermi momenta k F and k F , the 
neutron and proton densities n n and n p , and the total baryon density n = n n + n p , the 
potential energy density is given by 
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The single particle potential energy for a given nucleon with isospin index r and with 
momentum k is 
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(67) 



where rikk is the (rr)-diagonal element of the occupation number matrix 

(y\a+ai\^) , 



(68) 



each label i, j denoting momentum, spin, and isospin. For Skyrme interactions, one gets 
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The p-h interaction is obtained by functional differentiation of the single particle potential 
energy or equivalently the double functional differentiation of the total potential energy, 
namely 
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(k\k 3 \Vph\k4k2 ) 



6 2 (V) 



6n k3kl Sn k 



(70) 



Hereafter, we will employ the standard notation for the participating momenta, namely 
\t\ = q + qi, k 2 = q 2 , k 3 = q 1; k 4 = q + q 2 , which indicates that q is the transferred 
momentum. The p-h interaction can be expressed as 



^(qi,q 2 ,q)= £ V {r ^ q 2 , q) P {S) 

(rr',5) 



(71) 



with the p-h spin projectors P^ = 1/2 and P^ = a ■ a' /2. The symbol (tt') indicates the 
four isospin combinations 



(pp x ;pp x ) and (nn x ;nn x ) with r' = r, 
(pp' 1 ;™' 1 ) and (nn -1 ; pp -1 ) with t' = —1 



(72) 
(73) 

Two other isospin combinations, namely (pn -1 ; pn' 1 ) and (np' 1 ; np' 1 ), are needed to 
describe the p-h interaction relevant for the charge exchange processes. These are, however, 
not independent, and may be related to interactions in the (pp -1 ; pp -1 ), (nn -1 ; nn' 1 ), and 
(pp' 1 ; nn' 1 ) by isospin considerations. Denoting the p-h interaction in the spin independent 
channels as f pp , f nn , and f np for the particle hole states (pp~ 1 ;pp' 1 ), (nn -1 ; nn -1 ), and 
(pp' 1 -^^ 1 ), respectively, we arrive at the relations [fnj 
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(75) 



The interaction between p-h proton states f pp is related to f nn evaluated at proton fraction 
(1 — x) by isospin symmetry. The spin dependent p-h interactions characterized by g nn , g pp , 
and g np are obtained by taking functional derivatives of the energy density with arbitrary 
spin excess and are given by 



gnn = t (x - 1) + ^3^ 7 (^3 - 1) + 7 [tl{Xl - l) - t 2 (l + X 2 )) q 2 

o 4 
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+ r [tixi + t 2 x 2 ] (k F (n) 2 + k F (p) 2 ) . (77) 
o 

As noted earlier, the quantity g pp is equal to g nn evaluated at proton fraction 1 — x. 
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In the schematic models of Ref. [[Il]], one begins with a parametric form for the energy 
density for spin symmetric, but arbitrary isospin asymmetric, matter. The single particle 
potential, obtained by functional differentiation of the potential energy density, is given by 
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(78) 



where the upper (lower) sign in =p is for neutrons (protons) and % ^ j. Taking functional 
derivatives of the single particle potential energy density, we arrive at the p-h interaction 
parameters: 
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Their explicit algebraic forms are: 
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The form of the energy density employed in these schematic models does not explicitly 
account for the explicit spin dependence of the nucleon-nucleon interaction. The potential 
energy density is independent of any spin excess indicating that g nn and g np are zero. 

The p-h interactions parameters discussed in this appendix are related to the Fermi-liquid 
parameters. In symmetric nuclear matter the appropriate relations are 
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where the normalization factor Nq = 2M*kF/n 2 is the density of states at the Fermi surface. 
For neutron matter Fq = Nof nn and Go = SwA^o, with Nq = M*kp/TT 2 . 



APPENDIX B: POLARIZATION FUNCTIONS 

The various polarization functions required to evaluate the Hartree and RPA response 
functions are presented in this appendix. The zero temperature polarization functions may 
be found in Ref. |34| and for finite temperatures in Ref. []3"5|| . Here, we collect the present 



extensions of these results to asymmetric matter, and, in particular, to unlike p-h excitations. 
For space like excitations, q 2 , < 0, they are given by 

[m n l(jb , = */^ E ?-W 2 ™ 2t> e , (84) 
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In the above, 



6 = F(E* p , E; +q )[5(q - (E p+q - E p )) + 5(q - (E p - E p+q )\ , (88) 
F(E* p , E* p+q ) = f 2 (E;)(l - h(E* p+q )) , (89) 



E* = y/\p\* + M?, E p = E* p + U. (90) 
The particle distribution functions fi(E) are given by the Fermi-Dirac distribution functions 

fi(F*) = (1 + exp[(R* - v^/kT})" 1 , (91) 
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where v is the effective chemical potential defined by 

Ui = ^ - Ui = Hi - (g uBi wo + h Bi g P B&) , 



(92) 



and the particle labels 2 and 4 correspond to the initial and final baryons. 

The angular integrals are performed by exploiting the delta functions. The three dimen- 
sional integrals can be reduced to the following one dimensional integrals: 
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The lower cut-off e_ arises due to kinematical restrictions and is given by 
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where 
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It is convenient to re-express the polarization functions as follows: 
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These integrals may be explicitly expressed in terms of the Polylogarithmic functions 
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which are defined to conform to the definitions of Lewin []50| . This Polylogarithm represen- 
tation is particularly useful and compact: 
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We note that the nonrelativistic structure function for neutral current scattering, Eq. (|3|) is, 
aside from the factor M| T/nq, equal to Jo since £1 = In the case of neutral currents, 
some of the terms in the above are simplify, 
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and we recover the results obtained earlier in Ref. [55]. 

For the real part, analytic forms do not exist. In the following we present, in integral 
form, the required polarization functions. The results presented below generalize the results 
of Ref. |35| to the case when the particle and hole states correspond to different baryons. 
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The particle labels % and j refer to the initial and final state baryon. Anti-particle contri- 
butions have been neglected in the above formulae as they provide negligible contributions 
for the temperatures and densities of interest. Since the typical temperatures never exceed 
a few tens of MeV's they will be only important for electrons and only when fi e <C T. 

The vacuum contributions, which only play a role in the relativistic Hartree approxima- 
tion, are not explicitly given here, but may be found in Ref. |§. 
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FIGURE CAPTIONS 



FIG. 1: The ratio of the neutrino mean free paths computed with (M*) and without 
(M) effective mass corrections as a function of nucleon density n# in units of the nuclear 
equilibrium density uq = 0.16 fm~ 3 . Results are shown for pure neutron matter at T = 5 
and 30 MeV. The dash-dot curve shows the density dependence of the neutron effective 
mass. 

FIG. 2: The density dependence of the Fermi-liquid parameters and the neutron effective 
mass for neutron matter for Skyrme models studied in this work. 

FIG. 3: The density dependence of Fermi-liquid parameters of microscopic calculations 



from Refs. [28| and p9l (see text for details). 



FIG. 4: The neutrino differential cross sections in neutron matter as a function of the 
ratio of energy to momentum transfer for T = MeV (left panels) and T = 10 MeV (right 
panels) at n (top panels) and n /4 (bottom panels). Results are for the Skyrme model 
SLy4 for q = E v = 25 MeV. The symbols V and A refer to the vector and axial vector 
contributions, respectively. Thin lines refer to the Hartree results and thick lines to the 
RPA results. 



FIG. 5: Same as Fig. 4, except that results are for the model of Backmann et al. |29| and 
are shown for three densities: 2n (top panels), n Q (center panels), and no/4 (lower panels). 

FIG. 6: The density dependence of the neutrino mean free paths at energy E v = nT 
in neutron matter for T = 10, 20 and 30 MeV. Results for the microscopic calculations of 
Backmann et al. are in the top panels and those for the Skyrme model SLy4 are in the 
bottom panels. The right panels show the extent to which RPA corrections modify the 
Hartree results. 

FIG. 7: The density dependence of the neutrino mean free paths in the relativistic models 
RHA and GM3 in the Hartree approximation (i.e., RPA correlations were ignored) at energy 
E v = ttT for T = 0, 10 and 20 MeV. For reference, the results for the case in which effective 
mass corrections are ignored are also shown as dashed lines. 

FIG. 8: The neutrino differential cross sections in neutron matter in the model RHA for 
q = E v = 25 MeV. Dotted lines show the Hartree results, while solid, dashed and dot-dashed 
curves show the results including RPA correlations with g' = 0, 0.3 and 0.6, respectively. 
The left panels have T = MeV and right panels have T = 10 MeV; the top panels are for 
4no and bottom panels are for hq. 

FIG. 9: Same as in Fig. ||, but for the model GM3. 

FIG. 10: The density and temperature dependences of neutrino mean free paths at energy 
E v = tcT in neutron matter. The top panels are for the model GM3; the bottom panels are 
for RHA. The left panels show the results in the Hartree approximation. In the right panels, 
effects due to spin dependent correlations introduced through the Migdal parameter g' are 
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shown. The different curves correspond to the same temperatures as in the left panels. 

FIG. 11: The neutrino differential cross sections in symmetric nuclear matter for q = 
E v = 30 MeV. The Fermi-liquid parameters employed are given in Eq. (pop. Results for 
the free gas, the Hartree approximation, and with RPA correlations are compared for T = 
and 10 MeV. 

FIG. 12: The temperature dependence of the neutrino scattering mean free path in sym- 
metric nuclear matter at density no for the Fermi-liquid parameters in Eq. The left 
panel shows results for thermal neutrinos (E u = irT) calculated in the Hartree approxima- 
tion, and the right panel shows the effect of RPA correlations. 

Fig. 13: The density dependence of the effect of RPA correlations on the neutral current 
neutrino mean free paths in the nonrelativistic SLn2 models with different compressibilities. 

FIG. 14: The density dependence of the neutrino scattering mean free path in symmetric 
nuclear matter for the EOS SL22 at T = 10, 20, and 30 MeV. The upper panel shows results 
for thermal neutrinos (E u = ttT) calculated in the Hartree approximation, and the lower 
panel shows the effect of RPA correlations. 

FIG. 15: The density dependence of the Fermi-liquid and effective nucleon mass in 
symmetric nuclear matter for Skyrme models used in this paper. 

FIG. 16: The neutrino scattering mean free paths in symmetric nuclear matter for the 
Skyrme models SGII (top panels) and SLy4 (bottom panels), at a neutrino energy E v = 
7rT. Left panels show the Hartree approximation and right panels show the effect of RPA 
correlations. 

FIG. 17: Differential scattering cross sections in (3— stable neutrino-free matter for the 
schematic model SL22. The proton fraction x p = 0.049 for n# = n (lower panels) and 
x p = 0.22 at ub = 3n (upper panels). The left panels show results for T = while right 
panels are for T = 10 MeV. E v = q = 30 MeV are assumed. 

FIG. 18: Differential scattering cross sections in (3— stable neutrino-free matter for the 
Skyrme models SLy4 (top panels) and SGII (bottom panels), for the indicated proton frac- 
tions and temperatures. Results for both the Hartree approximation (dashed lines) and with 
RPA correlations (solid lines) are displayed. 

FIG. 19: The density dependence of the scattering mean free paths in (3— stable asym- 
metric matter for the schematic model SL22 at the indicated temperatures, in the Hartree 
approximation (left panel) and with RPA correlations included (right panel). In the right 
panel, thick lines show results obtained by employing empirical values for the spin parame- 
ters. Thin lines correspond to results obtained without any spin-dependent p-h interaction. 

FIG. 20: Differential cross sections at = 0.32 fm -3 in f3— stable neutrino-free matter 
for the field-theoretical models RHA (left panels) and GM3 (right panels) for T = (lower 
panels) and T = 10 MeV (upper panels). Results for the Hartree approximation alone and 
also including RPA corrections with g' = 0, 0.1 and 0.6 are compared. 
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FIG. 21: The density and temperature dependences of the neutral current mean free 
paths for (3— stable neutrino free matter in the field-theoretical model GM3. The upper left 
panel shows the Hartree results for the case E v = 3T. The influence of the spin correlations 
introduced via the Migdal parameter g' is strong, as can be deduced from the results shown 
in the upper right and bottom panels. 

FIG. 22: Charged current diferential cross sections in symmetric nuclear matter at nu- 
clear saturation density, for the case E v = q = 50 MeV. Empirical values of the Fermi-liquid 
parameters are employed, and results for the Hartree approximation and with RPA correla- 
tions are compared. For T = (left panel), the peaks labelled 'GD' and 'GT' correspond to 
the giant dipole and Gamow- Teller resonances, respectively. For T = 10 MeV (right panel), 
the collective states are significantly broadened and overlap. 

FIG. 23: Charged current differential cross sections in (3— stable matter for Yl = 0.4 and 
T = 10 MeV at — no (left panel) and % = 3n (right panel). The schematic model 
SL22 was employed. Results for the Hartree approximation and with RPA correlations are 
compared for the case q = E v = //„. 

FIG. 24: The density and temperature dependences of the charged current neutrino 
mean free path in (3— stable matter for the SL22 model assuming Yl = 0.4. Results for 
the Hartree approximation (left panel) are compared with those including RPA corrections 
(right panel) with g' = 0.6. 

FIG. 25: Same as Fig. 21, except for neutrino-free matter and for q = E v + \x e and 
E v = irT. 

FIG. 26: Same as Fig. 22, except for neutrino-free matter and for q = E v + /i e and 
E v = irT. 

FIG. 27: Charged current differential cross sections in (3— stable matter for Yl = 0.4 and 
T = 10 MeV at n B = n (left panel) and n B = 3n (right panel). The relativistic model 
GM3 was employed. Results for the Hartree approximation and with RPA correlations are 
compared for the case q = E v = //„. 

FIG. 28: The density and temperature dependences of the charged current neutrino 
mean free path in (3— stable matter for the GM3 model assuming Y L = 0.4. Results for 
the Hartree approximation (left panel) are compared with those including RPA corrections 
(right panel) with g' = 0.6. 

FIG. 29: Same as Fig. 27, except for neutrino-free matter and for q = E v + fx e and 
E v = nT. 

FIG. 30: Same as Fig. 28, except for neutrino-free matter and for q = E v + /x e and 
E v = tcT. 

FIG. 31: RPA corrections to the diffusion coefficients in the density-temperature plane 
for the field-theoretical model GM3. Results are shown for the three diffusion coefficients 
D 2 , Z?3, and D± (upper, middle, and lower panels, respectively) for neutrino- free matter (left 
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panels) and neutrino- rich matter (right panels). 



FIG: 32: The density and temerature dependence of the ratio of the RPA and Hartree 
diffusion coefficients (Df- PA / 'Df) in neutrino-free matter. 

FIG: 33: Central values of thermodynamic quantities (entropy s, temperature T, baryon 
density tib, v e concentration Y u , v e chemical potential fj, u , electron concentration Y e ) during 
the evolution of a baryon mass 1.4 M protoneutron star. The relativistic model GM3 for 
nucleons-only matter was employed. Results are compared for the Hartree approximation 
(baseline case) , the inclusion of RPA correlations, and a constant factor of three reduction 
of the opacities for densities greater than u = ns/n^ = 1/2 and 1/4. The latter case 



corresponds to that studied in Ref. | l2fl . 



FIG 34: The upper panel shows the total emitted neutrino luminosity for the protoneu- 
tron star evolutions described in Fig. 33. The curves are labelled as in Fig. 32. The 
lower panel shows the ratio of the luminosities obtained in the three models which contain 
corrections to the baseline (Hartree approximation) model. 

FIG 35: The total emitted neutrino luminosity for long-term protoneutron star evolutions 
described in Fig. 34. The curves are labelled as in Fig. 34. 
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FIG. 1. 
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FIG. 6. 
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